Extreme value statistics in coupled lasers 
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2^ Experimental configuration for investigating the dynamics and the statistics 

of the phase locking level of coupled lasers that have no common frequency 
is presented. The results reveal that the probabihty distribution of the phase 
locking level of such coupled lasers fits a Gumbel distribution that describes 
the extreme value statistic of Gaussian processes. A simple model, based 
on the spectral response of the coupled lasers, is also described, and the 
t/3 calculated results are in good agreement with the experimental results. @ 
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^ Phase locking of coupled oscillators was studied over the years in many different contexts 

including chemical oscillators with mutual coherence [1], arrays of Josephson junctions that 
are frequency locked [2] and arrays of coupled lasers that are phase locked [3-5]. In these, 
^ complete phase locking occurs when all the oscillators have at least one common frequency. 

^ When there is no common frequency, the oscillators group in several clusters, where each 

^ cluster oscillates at a different frequency [6]. So far there is very little, if any, experimental 

^ investigations which deal with the dynamics and statistics of coupled oscillators that do not 

have a common frequency. 
^ In this letter, we deal with the phase locking level in an array of lasers that have no common 

frequency, and show that the distribution of the phase locking level is in good agreement with 
the Gumbel distribution function which describes the extreme value statistics from Gaussian 
processes [7,8]. Specifically, we investigate the phase locking level of an array of 25 coupled 
fiber lasers. Although each fiber laser support 100,000 eigenfrequencics, the probability to 
find a common frequency for all the lasers in the array is very small (< 10~^) [9-11], so the 
lasers group in several clusters, each with its own frequency [12, 13]. Due to thermal and 
acoustic fluctuations, the length of each fiber laser and its corresponding eigenfrequencics 
changes rapidly and randomly. As we show below, phase locking minimizes loss in the array. 
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so mode competition will favor frequencies that maximize the size of the phase locked clusters 
at each moment [14] . Since the distribution of the phase locking level for different frequencies 
is Gaussian, the statistics of the maximum phase locking level is described by the Gumbel 
distribution function. 

The experimental configuration that we used for measuring the phase locking level of an 
array of fiber lasers is presented in Fig. 1 and described in detail in [12]. Briefly, each fiber 
laser was comprised of a Ytterbium doped fiber, a rear high reflection (> 99%) fiber Bragg 
grating (FBG), and a front low reflecting (5%) FBG both with a lOnm bandwidth. Each 
laser was pumped through the rear FBG with a 975nm diode laser at 200mW, and after 
the front FBG we attached a collimator to obtain a 0.4mm diameter beam. The collimators 
of all the 25 lasers were accurately aligned in a 5X5 square array of parallel beams with 
parallelism better than O.lmrad. The separation between adjacent beams was 3.6mm. A 
representative near-field intensity distribution, measured close to the output coupler when 
all 25 fiber lasers are operating, is presented at the lower inset of Fig. 1. We determined 
the length of each fiber laser by measuring the longitudinal mode beating frequency at their 
output by means of a fast photodetector connected to a RF spectrum analyzer. We found 
that the distribution of the lasers lengths is Gaussian with a mean value of 3m and a width 
of 0.5m. The intensity of each fiber laser was about lOOmW which is much above threshold 
but still low enough so nonlinear effects where negligible [15]. 



Fig. 1. Experimental configuration for phase locking an array of fiber lasers and for deter- 
mining their phase locking level. OC - output coupler, PR - partial refiector, BS - 50% beam 
splitter. Inset - near field intensity distribution when all 25 fiber lasers are operating. 

The coupling between the fiber lasers was achieved by means of four coupling mirrors 
denoted as ri, r2, and with refiectivity of 40% for ri and and refiectivity of 100% 
for r2 and r^. All the coupling mirrors were located close to the focal plane of a focusing 
lens with 500mm focal length, forming a self imaging cavity with the array. Since there was 
only enough space for one pair of mirrors within the Rayleigh range of the focussing lens, 
we inserted a 50% beam splitter to obtain another focal plane where we placed another pair 
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of mirrors. By controlling the orientations of the coupling mirrors we could realize a variety 
of connectivities for the fiber lasers in the array, and in our experiments we concentrated 
on the two-dimensional nearest neighbors connectivity [12]. Finally, we directed about 10% 
of the light with a partially reflecting mirror (PR), towards an output coupler (OC) of 99% 
reflectivity. The OC was placed at a distance of 2f from the collimator array, and reflected 
part of the light from each laser back onto itself with the same delay as the light that is 
coupled from the other lasers [16]. 

We measured the phase locking level as a function of time for different number of lasers 
in the array. This was done by continuously detecting the far-field intensity distribution of 
the interference pattern of all the light from the array with a CCD camera, determining the 
maxima and minima intensities, and calculating the average fringe visibility along the x and 
y directions. The fringe visibility provides a direct measure for the phase locking level that 
ranges from to 1. The correlation time of the phase locking level is shorter than 100msec, so 
over a 10 hours period we acquired about 370, 000 uncorrelated measurements of the fringe 
visibility. Representative experimental results of the fringe visibility as a function of time for 
10 seconds interval are presented in Fig. 2. The insets show two typical far field intensity 
distributions - one with low fringe visibility and the other with high fringe visibility. 
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Fig. 2. Typical experimental results of the phase locking level as a function of time over 
a 10 seconds interval. The phase locking level was determined from the far field intensity 
distribution of the output. Insets show typical far field intensity distributions - (a) low fringe 
visibility where the phase locking level is low; (b) high fringe visibility were the phase locking 
level is high. 

Using the 370,000 measurements, we determined the probability distribution of the fringe 
visibility and fitted it to the distribution of extreme values in Gaussian processes, namely 
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the BHP distribution [17]. A simplified form of BHP distribution is given by 




where a denote the mean value, /3 the width and c the measure for correlations in the Gaus- 
sian process, with c = 1.58 indicating highly correlated process and c = 1 an uncorrelated 
process. After fitting our probability distribution to the BHP distribution, we determined 
that c = 1.03, indicating that the Gaussian process in our case is uncorrelated. Consequently, 
the BHP distribution of Eq. (1), reduces to the generalized Gumbel distribution [7]. Repre- 
sentative experimental probability distributions for an array of 25 fiber lasers and an array 
of 12 fiber lasers with fits to a generalized Gumbel distribution functions are presented in 
Fig. 3. The inset shows the corresponding results in a linear scale. As evident, there is a very 
good agreement between the experimental probability distributions of the phase locking level 
and the Gumbel distribution function for the array of 25 fiber lasers, but not as good for the 
array of 12 fiber lasers. 
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Fig. 3. Probability distributions of the phase locking level for two arrays of fiber lasers. 
Asterisks (blue) - experimental probability distribution of the measured phase locking level 
of an array with 25 fiber lasers; circles (red) - experimental probability distribution of the 
measured phase locking level of an array with 12 fiber lasers; solid and dashed curves - 
Gumbel distributions fitted for both cases. Inset shows the probability distributions with 
the Gumbel fits in linear scale. 

First, we explain qualitatively the connection between the phase locking level and the 
Gumbel distribution. The number of lasers in each phase locked cluster changes rapidly and 
randomly, but always while maximizing the phase locking level in the array. This maximum 
phase locking level occurs at a specific frequency out of all the possible frequencies within the 
FBG bandwidth where the lasers losses are minimal. By considering the spectral response of 
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the coupled lasers it was shown, that the distribution of the phase locking levels for all the 
available frequencies is an uncorrelated Gaussian distribution [9,18]. Since the distributions 
of maxima of uncorrelated Gaussian processes are described by the Gumbel distribution 
functions [7, 8] , the probabihty distribution of the phase locking level should be the same. 

As evident from the results of Fig. 3, for the array of 25 fiber lasers the average phase 
locking level is 0.28. Accordingly, rare events with more than triple the average phase locking 
level can occur, so there is good agreement between the experimental probability distribution 
and the Gumbel distribution for phase locking levels up to 0.85. This is not the case for 
smaller arrays of fiber lasers where the average phase locking level is higher and the effect 
of clipping of the distribution at unity phase locking level is significant. To illustrate this 
effect, we measured the phase locking level for smaller arrays, with 20, 16 and 12 fiber lasers, 
and found that as the number of fiber lasers in the array decreases, the fit for a Gumbel 
distribution is less exact. This is clearly evident for the array of 12 fiber lasers where the 
average phase locking level is 0.57 and the agreement to the Gumbel distribution is only 
good for phase locking level below 0.7. 

Next, we present a simple quantitative model that relates the phase locking level to an 
extreme value in the spectral response of the array of coupled lasers. We start by assuming 
no gain and determining the spectral response of each laser cavity when we replaced all com- 
ponents that couple light into the laser cavity with an effective front mirror. The reflectivity 
of this effective mirror depends on the frequency [18]. For example, the effective reflectivity 
of the i'th laser in a one-dimensional array of N coupled lasers, when each of which is coupled 
to its two nearest neighbors, is 

1 - r(l - 2k) + Yz^^^e^M-^ + yz^^^mt^ 

where k denote the coupling to the two neighbors, r the reflectivity of the output coupler, 
li the length of the i'th laser, k the propagation vector, and and R^''^ the effective 
reflectivities from all the lasers above and below the i'th laser, given by 

R^^ = ^ -2 , (3) 

and 

^l* = ^ , ■ (4) 

and R^i"^ = R^^'' = 0. Now, introducing the laser gain together with a mode competition 
results in amplifications only at frequencies with high effective rcficctivity. 

We solved Eqs. (2)- (4) explicitly for an array of 25 fiber lasers and obtained the effective 
reflectivity of each laser as a function of the frequency. By counting the number of lasers 
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with effective reflectivity above a certain tliresliold, we obtained tlie number of lasers in the 
main cluster as a function of frequency which provided a direct measure for the phase locking 
level of the array as a function of frequency [17]. Representative results are shown in Fig. 4. 
It shows the probability distribution of the size of the main cluster together with a Gaussian 
fit for the tail of the distribution. The inset shows the size of the main cluster as a function 
of the frequency. As evident, the calculated probability distribution above the mean size of 
the main cluster is Gaussian. Then we selected the maximum size of the main cluster within 
the bandwidth of the FBG (lOnm) and determined the resulting phase locking level as the 
ratio of the size of the main cluster over the size of the array. 
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Fig. 4. Calculated probability distribution of the size of the main cluster together with a 
Gaussian fit for clusters larger than the mean size. Inset shows representative calculated 
results of the main cluster size as a function of the laser frequency. 

We repeated these calculations for 50,000 different arrays each with a different random 
fibers lengths. The fiber length of the i'th laser in each array was chosen to be + AZj, where 
li is the measured length of the z'th fiber and AZj a random length taken from a normal 
distribution with mean and lO/xm width. The results are presented in Fig. 5. It shows the 
probability distribution of the calculated phase locking level together with a fit to a Gumbel 
distribution. As evident, there is a very good agreement between the distribution of the 
calculated results and the Gumbel distribution, indicating that the effective reflectivity is 
suitable for modeling arrays of coupled lasers, and that the underlaying Gaussian process is 
the number of lasers in the main cluster. 

To conclude, we measured the probability distributions of the phase locking level for arrays 
of fiber lasers, and showed that they fit to Gumbel distributions. We also presented that the 
underlying Gaussian process is the number of lasers that have a common frequency. Our 
results can now be exploited to predict for arrays with an arbitrary number of fiber lasers 
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Fig. 5. Calculated probability distribution of the phase locking level from 50,000 different 
arrays of fiber lasers as well as a fit to a Gumbel distribution. 

what will be the probability to obtain a specific phase locking level. Finally, Finally, by 
operating the fiber lasers close to threshold we observe strong fluctuations also in the total 
power of the array which can be related to the statistics of extreme eigen values of random 
matrices [19]. 

This research was supported by the Israeli Ministry of Science and Technology and by the 
USA-Israel Binational Science Foundation. 
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